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ANTENNAS AND PROPAGATION 


UDC 621.396.96' 06 
ALLOWANCE FOR ANTENNA ERROR IN DETERMINING RADAR TARGET SCATTERING MATRIX 


Kiev IZVESTIYA VYSSHIKH UCHEBNYKH ZAVEDENIY: RADIOELEKTRONIKA in Russian 
Vol 26 No 3, Mar 83 (manuscript received after revision 10 Mar 82) pp 74-75 


[Article by B.A. Atayants and V.V. Yezerskiy] 


[Text] The target scattering matrix [1] 


< Si Sys 
== , | (1) 
Sn Sm 
in any polarization basis 
am e(P ot) ame(ts +712, - ‘ob? (2) 


where vob and T_, -- basis angle of orientation and ellipticity, 


pb 
respectively -- is done with antennas having two orthogonal polarizations 
which coincide with the unit vectors of the basis (2). However, because 

of the finite accuracy with which antennas are fabricated their actual 
polarization differs from the required figure. It is important to investigate 
the connection between discrepancies in measurements and discrepancies in 

the polarization parameters of the antennas, and to find a method to allow 

for errors. 


The voltage at the input of the receiving antenna can be found from the 
formula [2] 


~ 


V = ASB. (3) 


where the designates transposition, and the column matrices A and B 
describe the receiving and transmitting antennas. The elements of matrices / 











A and B are the projections of the radio wave formed by the antennas in the 
radiating mode onto the unit vectors e1 and e, of polarization basis (2). 


In the general case the errors in the antenna orientation and ellipticity 
angles differ for different polarizations. We shall designate the actual 
antenna polarizations (Wr» TA)» (Wes TO)» (Was t,) and ( B? tT)» respectively, 


for the receiving and transmitting antennas. Then 


AYA = VA— WE Sta = Uy bWa= Va Hp A= ttt, 
405 = ¥a— Yb» Sta =ta—T pbs Va =Ve— Hp, AtR= tat to,~ — (4) 


Knowing the results of measuring the voltage at the output of the receiving 
antenna Vane Vio» Vo9 for the three different antenna polarization 


combinations, we can obtain from (3) the followjng system of equations with 
respect to the unknown Sia? Sio S59 (assuming $4955): 


*e%,? 


V sq yb ,S,, + (apg Ogb,) Sy H+ gbyS 54, Vyq = Oyb/Syy + (03 “+ 36)) Syq + Ayb2S,0," 


- aia. tet. tele (5) 
Vig = 04b,S,, + (jb, + 49b;) S,, + agb75,,. 
The coefficients of Sip Si and Soo of this system can be found by 
using the unitary matrix obtained in [3] to calculate the transformations 
of the elements of column matrices A and B: 
a, = cos Arp’, cos At, + j sin Awp, sin (tat At,), a, = sin Aw’, cos ryt At’,) — 
— {cos AY, sin At, a, = — sin Ap, cos (At, — tH +- jcos Ay, sin At. (6) 


a; == COS AW, cos At, + j sin Ay, sin (At, _- “bp! 


bi» ba» bi» b, are described by analogy with (6), substituting the 


subscript B for A. 


By solving (5) for the unknowns § 9» we obtain the following 


11’ *42° 52 


in general form: 











San XanM a + Yan tg + Zan no 


te ; (7) 
hy mm 1,2, Xy me (— 1) tb) bg nl (d,by — 6765) (a,b, — @56;)], 
where 
Yi. an (— 1)*+*-! {b5_ nbs —_/((0,05 — byb5) (4,65 — a,b5)1 + 43 fyi, - a,a.,) x 
x (a,b,— a9b,)1}, 24, =(— 1y* +a, _4a3_,/(a,49 — Gay) (aib2 — @yb4)]. (8) 


These expressions are not valid when the determinant of system of equations 
(5) is zero, which can occur in one of these three situations: 


a 


aja,—aya,=0, bb, bib = 0, a’bs — ab = 0. (9) 


It is easy to demonstrate that one of these conditions can be satisfied 


in accordance with one of these practical cases: 1) Th at" at45°, i.e., 


both polarizations of the receiving antenna are circular and coincide; 

2) That! pats”, i.e., both polarizations of the transmitting antenna are 
° ' 7"! ° " ' ° 

circular and coincide; 3)T Kip ts or T A T pnts » 1eOoe Vi2 or Voy is 

measured with the receiving and transmitting antennas having coinciding 

circular polarizations. 


It is clear that these situations contradict the original statement of 

the problem of measuring the scattering matrix by antennas with different 
polarizations, and indicate that two of the three equations in system (5) 
become linearly dependent. Therefore, it is impossible for conditions (9) 
to be satisfied in practice. 


Expressions(7) and (8) thus make it possible to,calculate the scattering 
matrix elements in any polarization basis (e,> e,) by the results of three 
measurements of the complex voltages by antennas with polarizations which 
are different from those required. 


Let us look at three cases which are of practical importance. 


5 Ar," , = Art! 9 Ay," ° ss « w 
1. Aw A Aw A? AT A AT A? Aw B Aw B? AT B At B? i.e., the antenna 


errors are the same for both polarizations and differ for the receiving 





and transmitting antennas. It then follows from (5) that 


ted, amd, imdb = 
where * designates complex conjugation. 
Ze Ay" =A" p» Ay" =A" 5» At’ ,=At',, At" ,=At",, i.e., the same antenna is 
used for receiving and transmitting. Then 

' = " = ' = " ' = = " = ' es = " 
3. Aw A Aw A Ay B At B? AT A AT A AT 3 AT B? i.e., the conditions 
of the first two cases are satisfied, which yields 
ab =a,—6, a,=b,=—a=—b. (12) 


Consequently, the linear relationship (7) holds between the measurements 


of Vine Vior Yoo 11’ 532% S22 in the 


polarization basis (e,, @5)- The coefficients which enter into expressions 


and the scattering matrix elements S 


(7) are in general functions of the eight complex numbers a's, a's a" ys a", 
b')» b'., bY» b", (8), and they are expressed through the antenna polarization 


parameter errors (6). 


Expressions (8) can be simplified in certain cases by reducing the number of 
unknown variables to four (10), (11), or even to two (12). 


Obviously, these findings allow the a priori known measurement antenna 

' ' hi " ' ' Ai " 
parameters A? T A? Vea? TA? y Br. B? Wy B? T B and the results of the 
three measurements o! Viv Vi2 and Voo to calculate the elements of the 


scattering matrix in any polarization basis. 


This analysis is valid for single-position radar location, i.e., when the 
equality $9755) is satisfied. Violation of this condition makes it 
necessary to make four measurements of Ya Vio Vo4 and V59 for the four 


combinations of polarizations of the receiving and transmitting antennas 
and, accoringly, to solve a system of four linear equations instead of (5). 











The final expression for the sought S 
than the three in (7). 


* will contain four terms, rather 
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COMMUNICATIONS 


UDC 621.391.8 


CROSS SPECTRAL-CORRELATION PROCESSING OF SIGNALS IN DIFFERENT ORTHOGONAL 
BASES 


Kiev IZVESTIYA VYSSHIKH UCHEBNYKH ZAVEDENIY: RADIOELEKTRONIKA in Russian 
- Vol 26,No 3, Mar 83 (manuscript received after revision 28 Apr 82) pp 3-7 


[Article by N.V. Bebikh and A.I. Denisov] 


[Text] This study investigates aspects of cross spectral- 
correlation transformations in different orthogonal bases. 
A method is examined for computing the correlation function 
by the method of spectral-correlation transformation and of 
computing the spectra and correlation functions of certain 
signals; the structural diagram of a spectral-correlation 
converter is presented as well. 


Systems of non-sinusoidal functions have recently been used effectively in 
modern radio electronics. These include systems of Walsh, Haar and other 
functions which are convenient for digital processing by computer, as well 

as systems of Laguerre, Hermite and Karunen-Loew functions for investigating 
transient attenuating processes, etc. The choice of a system of basis 
functions is determined by the mean square approximation error, the amount of 
computational labor involved, the noise tolerance of the system of basis 
functions and the ease of hardware implementation. In order to provide an 
operational analysis of signals under study in the most suitable orthogonal 
basis it is helpful to have the capability to switch quickly between different 
systems of basis functions. 


The connection between spectral coefficients in two different systems of 
functions is demonstrated in general form in [1! 


Let us represent the investigated signal by a Walsh series and, define the 
Walsh-Fourier cross-transition matrix (M1: (M}=1/N[S, }x{w} , where 


(S.J -- matrix oi partial sums of trigonometric functions; (wt -- transposed 
Walsh matrix, with dimension NxN, N=2", where n -- order of system of Walsh 
functions. 


The elements of matrix [S, ] are determined by the following formulas: for 
the odd matrix rows 











2 mm/N 
Shm = On sin kxdx, 
2n(m—1)/N 


where k=1,2,..., N/2 -=- number of odd row in matrix; for even matrix rows 


, Qnm/N 
Sim = = cos Ixdx,' 
2n(m—1)/N 
where k=1,2,..., N/2 -- number of odd row of matrix; 2=1,2,..., N/2 -- number 
of even row of matrix; m=1,2,..., N -- number of matrix column; N=2" -- 


number of digitization intervals. 


The Fourier spectrum of the signal in question can then be computed by the 
formula 


N 
‘ 
i=l 
j=! 


where any -- spectral Fourier coefficient; a spectral Walsh coefficient; 
Mey -- Walsh-Fourier cross-transition coefficients; or, in matrix form, [A.J= 


= [M Jx{A.], where [A.J -- column vector of spectral Fourier coefficient; 
[A] -- column vector of spectral Walsh coefficient; [MJ -- Walsh-Fourier 
cross-transition matrix. 

The number of operations required to make the cross transition between 


spectral Walsh and Fourier coefficients can be defined by the recursive 
formula 


Lor =1+ 2" (2"'—}), 








where n=1,2,... -- order of system of Walsh functions; N=2" -- number of 
spectral coefficients. 


The Walsh-Haar cross-transition matrix (Mo and Harr-Fourier matrix [Moe 
can be calculated analogously: 


[Mon] = WN (SH) X (W)", (Marl = VN (Se) XX (HY, 


where (S,,] -- matrix of partial sums of Harr functions; (u}! ~- transpose 


of normalized Haar functions matrix, with dimension NxN. 


The number of operations required for the Walsh-Haar cross transition can 
be calculated by the formula 


Lwy = 1+ > "2" — 1). 


The number of operations for the Haar-Fourier cross transition is defined 
as Lape? (L724 (28 -1)(2™""-1), n=2,3..., for n=l, Luptl: 


The total number of operations required to compute the spectral coefficients 
in Walsh, Fourier and Haar bases with the help of the cross-transition 


matrices [Me] and [MJ will be: Lot =L,+Lyt+L,tL,, where L,=2" -= 


computation of partial sums of signal Yn question; L,=n2" -- calculation of 


Walsh coefficients with help of fast Walsh transform; L4=Lip -- calculation 


of Fourier coefficients with help of Walsh-Fourier cross transition matrix; 


Leb -- calculation of Haar coefficients with help of Walsh-Haar 


transition matrix. For N=16 Lot =16+64+71+156=307. 


W 
When a system of Haar functions is used as the initial system, the total 
number of calculations will be L =L_+L.+L,+L,, where L,=2" -- calculation 
tot, 1 ~2 3 °4 1 


of partial sums of signal in question; L =2(2"-1, -- calculation of Haar 


2 


spectrum with help of fast Haar transform; L -- calculation of Walsh 


3 HW 
coefficients with help of matrix [M5 L,=Lup -- calculation of Fourier 
coefficients by matrix. For N=16 L ot =16+30+71+329=446. 

H 


If the spectral coefficients are calculated in Walsh, Fourier or Haar 
bases without using cross transition matrices and with the help of fast 











transforms, the total number of operations will be Lote th ath ytl, » where 
4 


L,=2" -- calculation of partial sums of signal in question; L,=n2" -- 


calculation of Fourier coefficients with help of FFT algorithm (number of 


complex multiplications followed by addition and subtraction) ; L,=n2" -- 


calculation of Walsh coefficients with help of fast Walsh transform (number 
of additions and subtractions) ; L,=2(2%=1) -- calculation of Haar 


coefficients with help of fast Haar transform (number of additions and 
subtractions). 


The use of cross transition matrices to calculate spectra in different 

bases reluces the amount of computation, with the greatest gain being 
achieved when a system of Walsh functions is employed as the initial basis 
system. When implementing the device in hardware, in order to determine 

the spectral coefficients in different orthogonal bases it is best (in terms 
of hardware costs) to determine the spectral coefficients in one basis 
(e.g., a Walsh function basis), and then to switch to other function bases 
with the help of cross-transition matrices, which will consist of weighted 
resistance units and adders in terms of hardware. 


Figure 1 shows the structural diagram of the spectrum converter. The 
spectral coefficients from the output of Walsh spectrum analyzer 1 are 
input to the inputs of the weight coefficient units -- Fourier 3, Haar 4 
and Laguerre 5, the outputs of which are connected to switch 7. A signal 
from control system 8 causes the switch to connect the proper weight 
coefficient unit to adder unit 9, which outputs th- spectral coefficients 
in the proper orth»gonal basis in parallel. Converter 10 is used to obtain 
the serial spectrum. 


This device permits spectral-correlation conversion of signals. In order 
to do this, synthesized signal distribution unit 2 and multiplier unit 6 
must also be used. The correlation function K(T) in this case is computed 
as follows: 


N 
K(t) = > S)Sj4m, 
j=! 


m=0 


where m=0,1,2,..., N-l; j=1,2,3,...,N; 8, 


digitization interval; t=m/N -- shift parameter. 


-- value of synthesized signal in 


Since Walsh functions are piecewise constant over the digitization interval 
1/N, it is sufficient to determine the values of the correlation function 








at m/N points, assuming the changes between those points to be linear [2]. 
Distribution unit 2 connects to the inputs of multipliers 6 the required 
values of the synthesized signal, which are known in each digitization 
interval, as a function of the value of the shift parameter tT. The order 

of the switching is assigned by converter control system 8. With this method 
of obtaining the correlation functions there is no need to provide a time 
offset T between the multiplied signals (like in ordinary correlators, where, 
e.g. frequency-dependent delay lines are used for this). In the present case 
the time offset is assigned by simply changing the combination of pairs of 
factors from the values of the synthesized signal known in each digitization 
interval. 


When implementing spectral-correlation transformations by computer, it is 
best to represent the expression for the correlation function as follows: 


T 


r N N 
Kut) =a | Ould ~ > i ha t) Wal (i, t +t) df, 
0 =! 


(=| 


where a, -- spectral coefficients of expansion of signal x(t) into Walsh 


series; a, -- spectral coefficients of expanding the signal y(t) into a 


j 


Walsh series. 

The internal integral represents the cross-correlation function of the 
actual Walsh functions K (tT), and is independent of the type of signals 
Leing investigated; therefore, it can be calculated in advance for 


different values of T. 


Then the expression for the correlation function can be obtained in this 
form: 


N N 
Ky (t)= Nay NV a,Kw (0). 


j=l i=l 


If we have the pre-« i: 1lated matrices [K (tT) ] for various values of T 


and the column vectors of the spectral coefficients of the signals in 
question [A,], [Ai], we can calculate the correlation function by the 


matrix method in accordance with this algorithm: 
1. Calculation of column vector of inner sum: [o, (t) J=(K(t,) IxfA, J. 


2. Formation of auxiliary matrix [o(tT)] 


10 











oO; (T,) Oy (t;) «+e Oy (t,) 
(a (x)} az} 1%) Fa (2) ++. Oy (Ts) 


D4 (Tm) Fy (Tm) «Fy (Tm) . 


where j=1,2,...,N -- number of matrix row; m=1,2,...,N -- number of matrix 
colu-n. 
in out 


S(t) Sit) S(t) 


Ow 
é 
or 



























































3. Determination of correlation function matrix 
[K (t)] =[o (t)] X [Aj]. 


In the present case the matrices [K (tT) ] are calculated directly during 


computation of the correlation function; the intermediate results of 
computing [K(t)] for different values of T are not stored in memory, but 


are used immediately to calculate the current values of the inner sums of 
[o,(t)], and each subsequent matrix [K (tT) ] is written in the location of 


the preceding one, which reduces the computer memory size. 


ll 








A program for calculating the correlation function by the spectral- 
correlation transformation method has been written in “Assembler-BASS" 
language and has been implemented on an "Elektronika NTs03-T" microcomputer 
for N=16. A standard matrix multiplication program and matrix row 
calculation subprogram were employed. 


The results of calculating the spectral-correlation characteristics of 
some functions in Walsh, Haar and Fourier bases are shown in the form of 
time diagrams in Fig. 2, where S(t) -- investigated signal; ays Ay» a, - 





Walsh, Haar and Fourier spectra, respectively; S'(t) -- synthesized signal; 
K(t) == correlation function, 


Analysis of the findings makes it possible to assess the effectiveness of 
using a particular function basis for signal processing (by minimum mean 
square error, number of non-zero expansion coefficients, correlation 
coefficient, etc.,)i.e., it permits adaptation to the signals being 
investigated. This increases the effectiveness of quick analysis of 
complex signals and of synthesizing optimal correcting and control effects 
in automatic information processing and control systems. 
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UDC 621.391.2 


POWER SPECTRUM COMPUTATION ALGORITHMS FOR PULSED RANDOM PROCESSES 


Kiev IZVESITYA VYSSHIKH UCHEBNYKH ZAVEDENIY: RADIOELEKTRONIKA in Russian 
Vol 26 No 3, Mar 83 (manuscript received after revision 14 Jun 82) pp 7-12 


[Article by T.G. Pletneva, G.N. Rozorinov and S.D. Eydel'man] 


[Text] This study solves the problem of computing the 
power spectrum of a pulsed random process with determinate 
timing intervals in finite closed form. A computer program 
is described for calculating the power spectra of signals 
controlled by a regular Markov chain. A decomposition 
method is proposed which can be used to compute the spectra 
of digital signals formed from packets of varying length. 


Two-level digital signals with an assigned timing frequency are employed 
widely in modern communications, hydro-acoustic and computer equipment [1]. 


When deciding which of the existing digital signals is most suitable for 
effective transmission through a given propagation medium, it is important 
to be able to find the power spectrum of the signal. Knowledge of the 
power spectrum makes it possible to determine the amount of influence 
between bits within symbols. 


The limited bandwidth of a real data transmission channel causes the 
occurrence of high frequency and low frequency signal distortions. These 
distortions can be reduced by selecting signals which have the desired 
spectrum shape [2]. 


Two-level digital signals, which are mappings of fixed-length binary 
sequences (block codes), have attracted a great deal of attention 
recently in digital data transmission tasks. Each symbol in the input 
binary sequence is transmitted by a pulse with duration T, or by a 
group of pulses (a packet"). In this case the sequence of states s 


of the digital signal is modeled well by a process controlled by a 

finite ergodic Markov chain [3]. Signals formed from "packets" of 
different lengths are also important for practical purposes. The study 

[3] solves the problem of computing the power spectra of signals formed 
from "packets" of the same length which are used, e.g., in digital magnetic 
recording technology. 
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The present study proposes a method for solving the problem of calculating 
the power spectrum when the "packets" forming the signal are of different 
lengths. 


Formulas for calculating power spectrum. The pulsed random process in 
question consists of "packets" with duration rT, and the forms of the 
"packets" are defined (e =(€)), Egyrrees Eo) The process is controlled 


by a finite ergodic Markov chajn with period d and transition probability 
matrix P. We shall designate A=(A,,A grees A, ) the probability vector 





. oe iin Pe the Psa eae distribution of. the Markov chain in question; 


A) = lim Po, pS=pd SA, S81,2,.065 PO-1 -- identity matrix; 

nn? 
; ~(s, V) 5 US. we 

A= (A 19454, j=1, a’ 944 -- Kronecker's symbol; m p>p” e Ae; YaQ Licces 
(Y)_. 5Y (9) “ tT  -iut 

d-1; m‘'’=A_ P'e Ae ; =Ae Ae ys mee Ae, ; a=fe dt; r -- number of pulses 
uw oO Sy u 0 

with duration T in "packet"; q -- number of "packets". 


The formula derived in [3] for the averaged power spectrum of a random 
pulsed signal with "packets" of the same length controlled by an 
arbitrary finite ergodic Markov chain can be rewritten in the form 


d—| 


P (@) = ae D>p> P'Re( | — (P* — Ag) eidrory! +. 








L.vel yd 
Tm 
+ > \" § (@ — Qnh/drT) A | tien 7 
drT ° I 
__ | (1) 
sin (dr/2 — yr — T — 
+ Fain Una oa A art [ COS (11 — v) oT y2| Muy, 
where [1-(P9-A,)e ort) -- the inverse matrix of [-(P°-A eo SMT, 


Formula (1) provides the solution to the problem of calculating the power 
spectrum in finite closed form. It is useful in that it allows standard 
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programs for calculating the inverse matrix to be employed in computing 
the power spectrum. 


The next section examines a numerical algorithm for calculating the spectrum 
for the case of a random pulsed signal controlled by a regular Markov chain, 


which is the most important case in practical terms. In this case the 
formula appears as 


4jaty1 
P(a)= | 5 ' Myy COS (w — v) wT + 


pve! 
+- > mit cos (sr +p —v) oT + 
pve! se) 
+ mi COS (1 — v) oT ane of « y ar . 
Peme—9etl—a tor Meir )} 
Vel 


h=—@ 


q 
4% 
We recall that the discrete component in (2) is lacking if >» Artin =0, VE. 


j= 


In particular, this is the case when the matrix P is doubly stochastic 


q 
and 2c = 0 Vie. 


Algorithm for calculating power spectra of signals controlled by regular 
liarkov chain. Figure 1 shows the flowchart of the numerical algorithm for 
computing the spectrum in accordance with formula (2). The following 
notation is employed: 








' ~ ~ (S 
= > as cos (p —v) oT; Q,= S yin cos (sr + pp — v) wT; 
ph. vel p.vel sel 


r 


Lose + Sale My 


vel 4 hw—eo 
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The following is essential in the use of this algorithm. By using 
corollary (4.1.5.) from [4], we can demonstrate that the maximum number 
of required iterations Sa needed to ensure the required accuracy A of the 


spectrum computations is determined by the formula 





Sy = N [In A/In (1 — 2e) + 1), (3) 


where N -- power of matrix of transitional probabilities P in which all 
elements are positive; € -- smallest element of matrix pN (actually, the 
number of iterations is always much smaller). 
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(continued) 











Key to Fig. 1 (continued) 


ts ts 
13° Q. +Q, foi all w/w, 


if not, increase S 31 -- output [qvantity] 
ts 
32 -- end 


24 -- Is it the case that 67545793 29 -- Q 


25 -- input Ses 
26 -=- output [quantity] 


After Sy is determined, actions are taken to find the final probability 


vector A, i.e., the probabilistic solution of the algebraic system AP=A 
is found (e.g., by the Gauss method with selection of the main element). 
We note that if P is a doubly stochastic matrix, then A=(1,1,..., 1)/q. 
S=l1 is then assigned and the second term Q, of formula (2) is computed. 
If Sy7Sel, Q, is stored and then at and Q, are calculated, followed by 
the spectrum P(w/w,.) for all given W/W,» 


If Syfl, an exhaustive search 1<S <5. is made, potipSp is first computed, 


and then the procedure described above is repeated cyclically. 


If the signal in question has a large transition probability matrix, Ses 
may be large. In most practical cases the quantity See can be bounded. 


Then the algorithm in question can be used to calculate the maximum error 


SH for a given Se in accordance with the formula 


- ts oP ts 
64 = 1(Q, - 9)/0 |, (4) 


where QF -~- value of Q, in previous iteration. If the error is not 


excessive, Sis is increased no further. 


In the below example of calculating the spectrum by the algorithm described 
above, the null members of the matrix P include those elements whose values 
do not exceed 10°*. For Seg the calculation took 1.33 minutes, and for 
Sg 700 -- 12.06 minutes. 


Decomposition method in computing power spectra of signals formed from 
"packets" of different lengths. The above referred to the case in which the 
"packets" from which the signal was formed were of the same length rT. A 
whoie group of signals which are promising for use in digital magnetic 
recording and which are obtained by means of group recoding of the input 
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binary information in which a particular number of zeros or ones may not 
occur in a row [5,6] do not have this property. These signals include the 
(0,1)-code or Gabor code [5], the (4,9)-code or Franaszek code, the M2 FM 
[6] and ChMP-1 signals (Fig. 2, a-d). 





















































































































































(2./)—code (49)—code M°EM ChMP=1 
s Flr La? g o 
e/ Sy Sia = Sy 7 S7- 
“ft 4] | 
nl SiS, TO PY 
Lg O\/ c x) 
Soy Soa uJ u Sy SA Syy S32 —_ 
1/1/70 / if /; 
Sy Ssok- A ae Sy 5, aa 
/ yaw Se S¢05u ia ww d 
SoSoo 7 << Z|? f ' 
S7/Syo {- f 50550553 J oft 
2 4 r J7, - c) 
a) b) 
Fig. 2 


The power spectra of these signals are successfully calculated by the 
decomposition method, which allows these signals to be put into 
correspondence with a certain Markov chain which controls the decomposed 
random process. The possibility of this decomposition follows from the 
fact that the causal connections which determine the structure of the 
process for the class of pulsed random processes with determinate timing 
intervals in question correspond precisely to the probabilistic connection 
between its elements. 


The "packets" from which the signals in question are formed are thus 
divided into a definite number of "subpackets" of equal length, the 
probability of occurrence of which as part of its own "packet" is unity. 
The dimensions of the transition probability matrices, of course, are 
greater. 


The application of the decomposition method is illustrated by the example 
of computing the averaged power spectrum of the signal whose structure is 
explained in Fig. 2b. 
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Example of calculating power spectrum. The (4,9)-code (Fig. 2b) is a 
promising one for use in high density digital magnetic recording. This 
signal has a rather complex structure. By using the decomposition me thod 
a (4,9)-code can be represented by 24 alternating "packets". The "packet" 
length is Tr r=6 (T,=6T), q=24, and the transition probability matrix 





The matrices e,-e P P. are shown in Fig. 3. 
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It is easy to note that a transition, e.g., from "subpacket" Si to 


"subpacket" S_,,has a probability of occurrence of unity, since both of 
them are part of the "packet" Si» while the probability of a transition 
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from 'subpacket" S19 (which is part of "packet" s,) to 'subpacket" So4 
(which is part of "packet" Ss.) is 1/8. We are easily convinced that the 


pulsed random process which is obtained is controlled by a regular Markov 


chain. In order to use formula (1) on a computer, the values of A and 4 / Wp 
must be added to the data above. We shall assign A=1073, and w/w 0-6 in 
steps of 0.1. 


The results of calculating the averaged power spectra are shown as plots 
in Fig. 4, where the values of P(w), multiplied by the coefficient n° /2T,, 


are laid off along the ordinate, and the values of the normalized frequency 
w/W,, are laid off along the abscissa. 
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Fig. 4 


Discussion of results. Conclusions. It is apparent from the plots in 
Fig. 4 that the maximum spectral energy of the (0,1)-code (curve a) lies 
near the normalized frequency w/w, =1, i.e., in the region of the center 








frequencies of the passband of a magnetic record-play back channel. The 
constant component in the spectrum of the signal is small. These 
characteristics of the spectrum allow the (0,1)-code to be recommended for 
use when the recording density does not exceed 256 bits/mm. 
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The (4,9)-code (Fig. 4b) has a very useful spectral shape from the viewpoint 
of increasing recording density. The maximum spectral energy of the (4,9)- 
code lies in the vicinities of the normalized frequency w/W.=0.25. However, 


the spectrum of the (4,9)-code contains a constant component twice as large 
as that in the previous case, which results in substantial distortions to 
the signal as it passes through the record-play back channel. 


The MPM signal (Fig. 4c) differs from the others in Fig. 4 in that there is 
no constant component at all. The maximum value of the spectrum of the 
M2FM signal occurs at the normalized frequency w/w,=0. 3. These characteristics 


of the M2 FM signal spectrum, which are helpful for increasing recording 
density, have resulted in its increasing use, especially for use with 
floppy disks. 


It should be noted that the spectrum of an M2 FM signal as calculated by the 
author's method and the spectrum of the same signal cited in [7] are 
exactly the same. 


Curve d in Fig. 4 reflects the spectrum of a ChMP-1 signal. The substantial 
level of the constant and extreme low-frequency components of the ChMP-1 
signal, which are not passed by the magnetic record-play back channel, 

cause significant distortions. The smallest time intervals of a ChMP-l 
signal which corresponds to data zeros are practically not reproducible. 
This makes it much harder to use the ChMP-1 and similar signals in high 
density digital magnetic recording. 
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[Article by A.I. Rog and A.A. Sirota] 


[Text] A method is presented for estimating the error 
probability in signal recognition against the background 
of correlated noise with incomplete description of the 
correlation connections of the features. The proposed 
method is used to substantiate the possibility of 
replacing optimal decorrelation transformation by 
partially decorrelative discrete spectral transformations. 


The problems which arise im a number of practical applications, such as 
automatic recognition of pulsed signals of various forms: S,(t)=m, (t)+N(t), 


i=1,L, where m, (t) -- determinate signal component; n(t) -- normally 


distributed noise with correlation function o°R(T), involve using decision 
rules of the form [1,2] 


k => min ((x — m)"M~' (x —m)] => x€ Og. (1) 


In (1) X=(S(ty) see, S(t, _1)) -- vector of amplitude (feature) samples of 

> 
signal S(t) obtained in the time interval [0,T]; m -- vector of mathematical 
expectancies of features; M -- covariation matrix of features, 


M=||a, 01], a. o=0°R(|k-2|T/N). 


The use of optimal decision rules with structure (1) involves definite 
difficulties which arise in estimating the matrix M and which also 

require a very great deal of time to implement the matrix transformations in 
(1) during the decision making stage. It becomes necessary to examine a 
special class of quasi-optimal recognition algorithms based on a certain 
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"disregard" of the correlation between features, which in the present 
case is manifested in that a matrix K which is slightly different from 
the actual matrix M is employed in (1). 


Let us now examine the losses which occur when this class of quasi-optimal 


algorithms is employed. If the error probability of pairwise separation of 
signals for optimal algorithm (1) is defined by the familiar formular [2] 


~ ° 
| 7 
— |e di, 








J = ~- ~~ — 
Pl = @ (a4), ay) = 0,5 } (m, — m,)' Mo" (m,—— my), (D(x) = "er 


«/ 
A 


we find this probability for the quasi-optimal algorithm from the 
distribution parameters of the logarithm of the likelihood ratio of the 
two signals [2] with matrix K substituted instead of matrix M: 


In (%)] =x" K7! (m, — my) — & (rm, + my)" K™ (m, — my). 


It can be shown that 


M, {{In 1 (x)\/0,} = + 0,5(m, —m)"K~'(m,—m):; k= i}: 


D{ {Inf (x)/0,} = (m, —m,)"K“"MK~ (m, —m), 


and the probability of error in pairwise separation in this case is 


== _ T = + -* 
(m, — m,) K—' (m, —m)) 





Pi; =D(6,), By = 





2 \/ (m, —m;)"K~"MK~ (m, — m,) . 


The relationship between the probabilities PY and i is determined by 
the relationship between the quantities Oy and "44 which can be represented 
>_> + 


in the form below by introducing the difference vector a=m,—Mm, the vector 


j 
bek7 1/23 and the matrix p=k71/2yx-1/2, 








/ = ~!J | 1 tie _ —- 
a, =0,5 V a’K ®K?M"K?K ta=0,5 FB", 





o -t 44 .f, Lt Lt a ——- 
By=a'K °K *al2 |/ a’K ?K ?MK ?K *a=6'b/2" 6’po. 
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In order to estimate the losses occurring with the use of —_ class of 


quasi-optimal algorithms we examine the quantity Ama, 4-8, ij a, ,U- -min(f 45/43 )] 


for a given error probability for optimal algorithm Pp? and corresponding 


value of a 


ij° 
Using Bekkenbakh and Kassel's inequality [3], we obtain 


ij 


~~ > _\ me OM . (Ao 4 Ay)? 
a?2_/B?, = (b" Bb) (67 B™'by(b"b)? << Bown’ (2) 





where roy -- minimum and maximum eigenvalues of the matrix 


B=K™ 28,1934 We obtain from (2) the following formula for the losses 
occurring with incomplete allowance for feature correlation in the 
quasi-optimal recognition algorithms: 








AP = max P!;—(a,,), max P= (22,,. Viiwo ) (3) 
ho + An—/° 


Thus, for an algorithm which does not allow fully for sample correlation 
(K=L, where I -- identity matrix) 


k => min {(x —m,)" (x— m)] => x€ 0, (4) 


we can write the following expression for the maximum possible recognition 
error probability (P y7const): 





2) ary) ~1/p0 
pi =0(a, : , a,,== D* (P3), 
j 1 Ke F Ryo ij ij (5) 
whe re Ko» Kye1 7 eigenvalues of matrix M. 
The difficulties involved in using information about the correlation of the 
initial features -- signal amplitude samples -- can be avoided by using 


the familiar decorrelation method. An optimal decorrelative transformation 
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is the expansion with respect to discrete eigenfunctions of the matrix M 
using a Karunen-Loew expansion [2]. The optimal decision rule in the 
space of the new uncorrelated features assigned by the vector y appears as 


k=> min |((y —n)'Q" (y—n)) => y€ O,. (6) 


where y=Vx, Qevmv', n,=Vm, , V -- optimal decorrelative transformation 


(ODT) matrix. The main shortcoming of the ODT is the difficulty of finding 
the expansion functions and the excessive amount of time required to 

execute the ODT. It is therefore better to try to replace the ODT with 

some other orthogonal transformation, which may have slightly poorer 
decorrelative properties but which allows the transformation procedure to 

be speeded up and simplified significantly. Such transformations include 
Fourier, Walsh and Haar discrete spectral transformations (DST), for which 
"fast" transformation methods have been worked out. The idea of the ODT-DST 
replacement is based on certain facts of the similarity of the structure of 
the bases of orthogonal functions of the optimal decorrelative transformation 
and of spectral transformations assigned on a discrete finite interval [4,5]. 
For Ne, for example, [5] notes that the Fourier transform approaches the 
Karunen-Loew transform asymptotically. 


The problem of estimating the effectiveness of the use of an “approximate” 
decorrelative transformation instead of the ODT is equivalent to the problem 
of estimating the losses which occur with an incomplete description of the 
correlation connections of the signal samples in quasi-optimal recognition 
algorithms. Actually, replacing the ODT with some DST assigned by the 
matrix H actually consists of using in place of optimal algorithm (6) an 
analogous structure of a quasi-optimal algorithm implemented in the space 
of the spectral coefficients of the DST -- z: 


hep min[(2—1)"E7! (@ — I) 2€ Ons (7) 


> + > T 
where z#Hx, fH E=diagP, P=HMH -- covariation matrix of spectral 


coefficients. 


4° 


Algorithm (7) doves not make complete allowance for the correlation between 
the spectral coefficients, and is quasi-optimal with respect to an algorithm 
with structure (1); however, it is realized in a spectral coefficient 

space in which the covariation matrix of these coefficients P is replaced 
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by diagonal dispersion matrix E. We can therefore estimate the losses which 
occur when quasi-optimal algorithm (7) is used instead of optimal algorithms 
(1) and (6) on the basis of expression (3), which is fairly general. The 
expression for the maximum possible Recognition error probability for 
algorithm (7) for a given value of Pay -- the error probability of algorithm 


(6) -- which determines the effectiveness of the ODT-DST substitution, looks 





like 
h2 2) Kory —_ 
7 0(a, re). oj =D" (PH), _™ 


where Xo» Xyiy 77 eigenvalues of matrix p7l/2pp-1/2 
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Formula (8) was used as the basis for investigating the effectiveness of 

using Fourier and Walsh DST in quasi-optimal signal recognition algorithms. 
Figure 1 shows the probabilities of pairwise signal recognition errors for 
different quasi-optimal algorithms as a function of the ratio of the noise 
correlation time to the signal duration. For comparison, the error probability 
Pt was calculated simultaneously for quasi-optimal algorithm (4), which does 


not make complete allowance for sample correlation. Figure 1 a, b and c show 


pki and pk2 as a function of the ratio of the noise correlation time to the 


ij ij 











signal duration -- T/T for different values of N and types of 


correlation function R(T) for a value of p? =0,01 (the dotted line in 


ij _ 
Fig. la corresponds to N=32, and the solid line to N=16 for R(t)=e a|t] , 


In Fig. 1 b and c for R(t)*e"|T| cos wt and R(t)=e"*" , respectively). 
It follows from these relationships (1 -- probability PkL 2F, 2W -- 


probability hi 


that the losses which occur with the use of quasi-optimal algorithms with 
structure (7) in Fourier and Walsh coefficient space are significantly 
smaller than for quasi-optimal algorithms (4) realized in the sample 
space. The decorrelative properties of the dst are particularly strong 
for small values of T/T, when the quantity Pty is practically the same 


computed for Fourier and Walsh transformations, respectively), 


as ave As follows, e.g., from Fig. la, by comparison with algorithm (6), 


based on the use of the ODT, signal recognition accuracy for noise 
described by an exponential correlation function is no more than 3-5% 
lower when quasi-optimal algorithms (7) with Fourier transformation are 
employed, and no more than 5-10% lower for Walsh transformation for values 
of T,/T<0.4, 
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[Text] Spectrum analysis devices can be simplified significantly by 
approximating the transformation kernel with piecewise-constant functions 
with a small number of levels [1]. It is particularly attractive to do the 
approximation with the help of a periodic two-level signal: this makes it 
possible to increase dynamic range significantly and to simplify the 
analyzer structure. The spectrum of such a periodic two-level signal must 
contain a first harmonic, after which the series of harmonics (2,3,...,p) 
is lacking, and the amplitudes of the higher harmonics (p+tl, p+2,...,) 

are arbitrary. 


The objective of the present study is to synthesize a periodic digital 
signal whose spectrum does not contain harmonics 2, 3, ..., p. We shall 
represent the two-level signal which has the required spectral structure 
in the form 


N 
yith= Sea(t—U— ty). £€(0.7), 
| 


where h(t)=f2? 6 (0,t;) 


0, for other T -- describes the shape of rectangular pulses 
with duration Ch the spacing of which in the periodic sequence of 


pulses y(t) with period T=Nt, is determined by the values of the coefficients 


i 


E, = (9, 1). (1) 
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We represent the complex amplitude of the ith harmonic as 


T N 
1 7 2 /sin(ni/N SN. in 
0 j=] 


For i<N the condition C,=0 is equivalent to the requirement 


N 
> ge /2AU1/N 9 (2) 


lea] 


The simplest signals which correspond to the solutions of (2) are groups 
of m pulses spaced apart by N/(mi) pulse positions, i.e., by 1/m part of 
the period of the ith harmonic component (Fig. 1). The requirement of 

(2) is satisfied for the signal shown in Fig. 1, since the left part of 
(2) can be viewed as the sum of m vectors on a complex plane (Fig. 2) with 
norm of unity and arguments 2tTn/m, n=1,2,..., m (which can be obtained by 
assuming 2=nN/(mi) in (2)). 





a 
‘ q Z 
N/m ~ uv ow 
Fig. l Fig. 2 


We note that the condition C,=0 is equivalent to having coefficients a, and 


by of zero in the expansion of y(t) into a real Fourier series, since 2C = 


=a,-jb, [2]. 
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The other solutions of (2) can be obtained by linear combination of the 
simple solutions; however, (1) must not be violated when this is done. 


Let U, be a group of m pulses corresponding to the solutions of equation (2) 


for a given i. In a periodic signal this group repeats every N positions. 
In the spectrum S, (2m£k) of such a periodic signal U, (t) the amplitude of 


the ith harmonic is zero: S, (21£,1)=0, where fy=1/T -- frequency of first 


harmonic of signal y(t). The convolution spectrum W(t) of such signals U 


’ 
1=2,..4.65P . 


W (t) = Ug (t) a Usa (te # UD (3) 


p 
is equal to the product of the spectra of these signals S_, (2Afgt) = [] S, (2fo%) 


i=? 


and, consequently, the amplitudes of the harmonic components i=2,...,p in 
it are zero. The convolution operation in (3) is the discrete convolution 
N 
f : ; 
of periodic pulse sequences of the type Ai) = S BWC 2), j=l... N. 
cd 


The signal W(t) is also periodic, with period T=Nt Obviously, the number 


i- 
of positions N within a period of this signal is a multiple of all factors 
of the type m2, 1=2,...,P-. 


In general terms, this is the algorithm for constructing the sought signal 
W(t). It should be noted that the presence in the convolution (3) of a 
certain group U,, as will be shown later, makes the amplitude of the 
harmonic numberedi equal to zero, as well as a series of higher harmonics. 
This makes it possible to reduce the number of groups in the right part of 
the convolution (3) and, consequently, to reduce significantly the number of 
positions N in the period of the signal W(t), which is extremely important 
for practical purposes. 


Let us now examine the spectral structure of a T-periodic sequence of groups 
of m pulses spaced t=T/ (im) apart. Finding the frequencies of the harmonics 


which are lacking in the spectrum of this sequence involves solving this 
equation: 


F (fy (X™-' 4x"? 4...4X+41)=0, where yaeJ2TET/ (im) 
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and corresponds to a delay operator by T/(im), while F(f) is determined 
by the shape of the pulses. 


Obviously, the equation 


XMV XM Xf Lm (X™ — IX — 1) = 0, (4) 


j2mk/m 101,2,... mel. 


There corresponds to each of the solutions 95 {A} a series of fregyerine: 
which are found from the formula X =X=ed "£T/ (im . However, xX =e ™k/m 


is a special case of (2) and has the solutions X =e 


consequently, 2mfT/(im)=27k/m+2mr, r=0, +1, +2,..., whence 


f=(k-+-rm)i/T, k=1,2..,m—1,7=0, 41,42, (5) 


We note that the factor (k+rm) in the right part of (5) can take on any 
integer values except for integer multiples of m. Consequently, the 
presence of a group U; consisting of m pulses in the convolution (3) leads 


to compensation in the resulting signal W(t) of all of the harmonics whose 
numbers are multiples of i, except for those whose frequencies are 
represented in the form fo=f im, £=0,1,2,.... In order to construct a 


signal W(t) in which the amplitudes 2,3,...,p of the harmonics are zero 
it is necessary to perform a convolution of the type (3) which includes 
groups U, of m pulses such that i runs through all of the primes and all 
of the integer powers of m not exceeding p. Let I(p,m) denote the set 
consisting of these values of i. As was noted earlier, the number of 
positions N within a period of the signal W(t) must be a multiple of all 
factors of the type im. Therefore, the minimum N is defined as N=mLCM 
{I(p,m)}, where LCM -- least common multiple. 


As m increases the number of groups U, in convolution (3) becomes smaller. 


i 
However, the resultant signal W(t) remains a two-level signal only for 
m=2 and m=3. 


The highest suppressed frequency f =foP is in principle unbounded. However, 


h 
for large p(p>10) the required number of positions N in a signal period 
increases immeasurably. 


Table 1 shows the minimum N=N(p) for the case in which the signal W (t) 
is synthesized by convoluting simple two- and three-pulse signals 
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(for m=m,.=2 and m=m 


2 373) and q -- number of pulses in period of W(t). 














TABLE 1 

P alma] N4/Ns | 92/%s | p Imm N,/Ng 92/s 

2 2/3 4/6 2/3 9 —/3 1680/1890 64/243 
3 2/3 12/18 4/9 10 —/— 1680/1890 64/243 
4 2/ - 24/18 8/9 1) = 2/3 18480/20790 128/729 
5 2/3 120/90 16/27 2 —/— 18480/20790 128/729 
i im 120/90 16/27 13 2/3 240240/270270 256/2187 
7 2/3 840/630 32/81 14 —/— — 240240/270270 256/2187 
Ss 2/— 1680/630 04/81 5 —/— 240240/270270 256 /2187 


lo 2/— 480480/270270 5912/2187 


The class of two-level signals W, Ce) can be used for spectrum analysis 


<f< 
of signals occupying a limited frequency band = . f — It can be 


i 
shown that there is no methodical analytical error if p>f_ /f_, -l, 
max min 


which confirms the advisability of adapting the synthesized signal W(t) 


to the analyzed process. One feature of this approach, in contrast to 
universal algorithms (specifically the FFT [3,4] where the transformation 
kernel is not approximated) is the fact that by increasing the number of 
discrete samples of the analyzed signal (increasing the digitization 
frequency) it is possible to eliminate the multiplication operation 
completely. In order to obtain the expansion coefficients of the 
analyzed signal it is sufficient to add the values of the samples taken 
at the instance corresponding to non-zero levels of the digital signal 
W(t). In calculating different expansion coefficients of the signal 


in question, the samples can be added independently and in parallel. 


Increased speeds can be expected in analyzing signals with relatively 

narrow bandwidths, especially when the code word lengths in the computer 

are long. It should also be pointed out that it is possible in principle 

to filter the analyzed signal first in order to represent it as the sum 

of narrowband components, and then to apply the analysis algorithm in question. 
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EFFECTIVENESS OF DIGITAL TRANSMISSION EMPLOYING MODIFIED WAGNER METHOD 
TO CORRECT SINGLE ERRORS 


Kiev IZVESTIYA VYSSHIKH UCHEBNYKH ZAVEDENIY: RADIOELEKTRONIKA in Russian 
Vol 26 No 3, Mar 83 (manuscript received after revision 23 Sep 82) pp 70-72 


[Article by R.E. Gut, M.Ya. Lesman and A.S. Shutov] 


[Text] Various decoding procedures which allow for the reliability of the 
decisions made about the signal elements transmitted are employed 
extensively in digital transmission systems. For transmission over radio 
channels these procedures are especially convenient when multi-channel 
modems are used (such as the Kineplex, MS [expansion not given], etc. [1]). 
In contrast to the approach usually taken [2], it is best to judge the 
symbol reliability according to the estimates of the signal/noise ratio 
provided by meters which are specially added to each of the modem channels. 


The present s.udy solves the problem of estimating the noise tolerance of 
parallel reception of signals coded with an (M, M-I) code and decoded by 
the Wagner method. The indicators which determine noise tolerance are the 
probabilities of correct reception of a code combination and a symbol in an 
arbitrary phased channel. Channel fading is assumed to be statistically 
uniform and independent. 


It can be demonstrated that the expressions for the mean probabilities of 
correct reception of a code combination Q, and of a symbol in an arbitrary 


jth channel Q, in the case in question appear as 


7 @ oo re) M—! 
=a YT omnis, h) fF uve F(S|h)dh\ as, = 
0 0 


where p-- mean error probability during reception of elementary signal 

in any of M statistical uniform channels; p(h) -- current error probability 
in receiving elementary signal as function of actual signal/noise ratio; 

M -- number of symbols in code word (number of parallel channels); f(h) -—- 
density function of signal/noise ratios; f(S,h) -- combined density 
function of estimate of signal/noise ratio and actual value of signal/noise 
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ratio; F(S|h) -- conditional distribution function of signal/noise ratio 
estimate assuming that the actual signal/noise ratio is h2; 


ad = M—1 
Q = (=) | 1s) {eH (s) —[ FS) —2 | per F (Sm) aa | dS -+ 
0 0 


(2) 


+ | FM-1 (5) { p(n) f(S, h) dhds, 
0 0 


where f(S) and F(S) -- unconditional density function and distribution 
function of estimate of current value of signal/noise ratio, respectively. 
















































































Fig. 1 Fig. 2 


Let us now specify the formulas derived above for the case of optimal 
incoherent reception of signals employing single phase-differential keying 
[2]. We shall assume that the channel fading is described by a Rayleigh 
distribution, and that the receiving location employs the devices for 
measuring the current signal/noise ratio which are described in [3], for 
which 


f(S\h)= ({1/T (v)] (4vh2)¥S¥—! exp (— 4vh?S), (3) 








where 2v -- size of sample used to form estimate of signal/noise ratio. 


By substituting (3) in (1) and (2), we obtain the following analytical 
formulas: 




















| (4) 
ea i . po 
2+) 2 J+ Rd—otl easmrmacer | | 
9 1 
l he -_ 
tee fem ft, 
2(1 + h2) 2 (1 Ad) (1 2 (1 — x)” or 
1 
A vy x¥M—Idx 
"2 [1 + Aj (lt — xt! 


Calculations were done in accordance with (4) and (5), some of the 
results of which are shown in Figs. 1 and 2. Figure 1 shows (1-Q,) 


during reception of a code combination as a function of the sample size 
(n=2v) for M=2 (dotted lines) and M=6 (solid line) and for two values of 
hy “=10 (curves 1 and 3) and h)*=100 (curves 2 and 4). Figure 2 shows 


the same for (1-Q,) with M=6 and hy =10 (curve 5) and hy =100 (curve 6). 
We note that when m=2 the probabilities Q, and Q, are the same. Furthermore, 


the straight dot-and-dash lines in Figures 1 and 2 show the values of the 
corresponding error probabilities during code combination and symbol 
reception which are achieved in a hypothetical system belonging to the class 
in question in which precise information is available about the actual 
signal/noise ratios in each channel. 


Analysis of the findings permits the following conclusions: 
1. The use of the modified Wagner procedure (considering the signal/noise 
ratio meters) allows the error probability for code combination and symbol 


reception in an arbitrary channel to be reduced by almost an order of 
magnitude, even in statistically uniform channels. This procedure should 
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provide even greater effectiveness when one of the channels is affected 
by spectrally concentrated noise. 


2. The mean error probabilities for code combination (1-Q,) and symbol 
(1-Q,) reception in an arbitrary channel are practically the same as the 


corresponding figures for the hypothetical system with a sample size of 
n=2 =15-20. 
V 
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